Collective Dipole Oscillation of a Spin-Orbit Coupled Bose-Einstein Condensate 
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In this letter we present an experimental study of the collective dipole oscillation of a spin-orbit 
coupled Bose-Einstein condensate in a harmonic trap. Dynamics of the center-of-mass dipole os- 
cillation is studied in a broad parameter region, as a function of spin-orbit coupling parameters 
as well as oscillation amplitude. Anharmonic properties beyond effective-mass approximation are 
revealed, such as amplitude-dependent frequency and finite oscillation frequency at place with diver- 
gent effective mass. These anharmonic behaviors agree quantitatively with variational wave-function 
calculations. Moreover, we experimentally demonstrate a unique feature of spin-orbit coupled sys- 
tem predicted by a sum-rule approach, stating that spin polarization susceptibility-a static physical 
quantity-can be measured via dynamics of dipole oscillation. The divergence of polarization sus- 
ceptibility is observed at the quantum phase transition that separates magnetic nonzero-momentum 
condensate from nonmagnetic zero-momentum phase. The good agreement between the experimen- 
tal and theoretical results provides a bench mark for recently developed theoretical approaches. 



Many interesting quantum phases can emerge in solid 
state materials when electrons are placed in a strong mag- 
netic field or possess strong spin-orbit (SO) coupling, 
such as the fractional quantum Hall effect [l| and the 
topological insulator Q. In cold atom systems, albeit 
neutral atoms have neither charges nor SO coupling, the 
recent exciting experimental progress demonstrates that 
artificial gauge potentials can be synthesized in labora- 
tory by laser control technique |3l4lfj|. Synthetic gauge 
potential is becoming a powerful tool for simulating real 
materials with cold atoms. Moreover, the system of SO 
coupled bosons does not have an analogy in conventional 
condensed matter systems, and can exhibit many novel 
phases such as striped superfluid phase [12|, [l3[ and 
half vortex phase (l4- 17|. 



Collective excitations play an important role in study- 
ing physical properties of trapped atomic Bose-Einstein 
condensates (BEC) and degenerate Fermi gases. Collec- 
tive dipole oscillation is a center-of-mass motion of all 
atoms. For a conventional condensate, the dipole os- 
cillation is trivial: the frequency is just the harmonic- 
trap frequency, independent of oscillation amplitude and 
interatomic interaction. This is known as Kohn theo- 
rem [UEiL For a SO coupled condensate, however, it 
was found 4[ that the dipole-oscillation frequency devi- 
ates from the trap frequency and the experimental data 
thereby can be explained by effective-mass approxima- 
tion. Recently, much theoretical effort has been taken 
to understand dynamics of a SO coupled BEC US-HI, 
and many predicted unconventional properties remain to 
be experimentally explored. In particular, the so-called 
sum-rule approach predicts [25] a unique feature of SO 
coupled condensate: spin polarization susceptibility-a 
static physical quantity-can be measured via dynamics 
of dipole oscillation. 



In this letter we experimentally study the collective 
dipole oscillation of a SO coupled 87 Rb BEC, occur- 
ring both in momentum and magnetization. The oscil- 
lation frequency is measured along various paths in the 
phase diagram, as a function of SO-coupling parameters 
and oscillation amplitude. Anharmonic properties be- 
yond effective-mass approximation are observed, includ- 
ing amplitude-dependent frequency and finite oscillation 
frequency at place with infinite effective mass. The ex- 
perimental data fit well with variational wave-function 
calculations. Moreover, following the proposal by the 
sum-rule approach [25|, we deduce the spin polarization 
susceptibility from the amplitude ratio between momen- 
tum and magnetization oscillations; the results are in 
good agreement with theoretical calculations. In partic- 
ular, the measured spin polarization susceptibility does 
diverge at the quantum phase transition that separates 
magnetic nonzero-momentum condensate from nonmag- 
netic zero-momentum phase. To the best of our knowl- 
edge, this is the first experimental measurement of a di- 
vergent spin susceptibility at magnetic phase transition 
in cold atom system. 

The experimental layout is sketched in Fig. [Ha). A 
87 Rb BEC of about 2.5 x 10 5 atoms is produced in 
the crossed optical dipole trap with wavelength 1070nm, 
beam waist 80/im, and trap frequency uj — 2ttx 
{45,45,55}Hz. A bias magnetic field Bbias is applied in 
the z direction. The BEC is illuminated by a pair of 
Raman lasers in the x-y plane with beam waist 240/xm, 
relative angle 8 — 105° and wavelength A = 803. 3nm. 
The Raman lasers couple the three internal states of the 
F = 1 manifold as shown in Fig. [Tfb) . By setting the 
quadratic Zeeman shift e = 3.37-E r with recoil energy 
E t = fc 2 /(2m) = 2-7T x 2.21kHz (fc r is the recoil momen- 
tum), we effectively suppress state \itif = 1) and prepare 
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FIG. 1: (a) Experimental layout. Field Bbias is along the z 
direction and the Raman lasers propagate in the x-y plane, 
(b) Raman coupling scheme within the F = 1 manifold, (c) 
Single-particle phase diagram in the Sl-S plane. The disper- 
sion spectrum £(k x ) has two (one) local minima in the blue 
(grey) regime. Experiments are along paths Pi, P2, and P3, 
with £ (k m ) in the insets, (d) Experimental time sequence, (e- 
f) Dipole oscillation with Q = 3.3-E r and <5 = 1E T along x(e) 
and z(f) direction. For comparison, oscillations without SO 
coupling are displayed as the grey dashed curves. 



a spin-1/2 system by regarding state \m-p = — 1) as | f) 
and |tof = 0) as | 1). This leads to a single-particle 
Hamiltonian as (h = 1) 



(fc x +fc r r 



Hn = 



2m 



n 
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(1) 



where k x is the quasi-momentum, is the strength of Ra- 
man coupling, and 8 is the two-photon Raman detuning. 
Diagonalization of Eq. (p} yields two eigenstates for each 
k x , in which the spin and the momentum of an atom are 
coupled. For small Q and 8, energy dispersion £(k x ) has 
two local minima, and bosons condense in one of the min- 
ima for our current experimental conditions. For large Q 
or (5, £{kx) has only one minimum. The phase boundary 
between the double- and the single-minimum region is 
displayed in Fig. (He), in which the perturbative effect of 
state \mp = 1) has been taken into account (as well as 
in theoretical calculations later). We shall systematically 
study the dipole oscillation along paths Pi,P2, and P3 
in Fig. (He). 

The time sequence for the experiment is shown in 
Fig- Did). After the BEC is prepared in the trap (t = 0), 
the Raman coupling is adiabatically ramped up from zero 
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FIG. 2: (a) Oscillation period T x along path Pi (S = IE T ). 
versus Q.. Black triangles are without SO coupling, and 
the blue line denotes the estimated trap frequency. Purple 
squares are with SO coupling for relatively small (below 0.2fc r ) 
oscillation amplitude, while green circles are for large (about 
0.6fc r ) amplitude. The inset is for T z versus Q.. (b) T% ver- 
sus oscillation amplitude with (fi = 4.8, 8 — 1.0)E r . Black 
squares with error bars are experimental data, (c) and (d), 
T x versus Q. along path P2 and path P3. In all plots, the red 
lines are the theory curves by solving Eq. (J2|. 



to Q in a time period t\ from 70 to 100 ms, and Bbias 
is slowly ramped, which adiabatically changes 8 from an 
initial value 5i > 20£" r to an intermediate value i5 m . At 
t\, detuning 8 is switched from <5 m to 8{ (black line in 
Fig. [Ijd)) in 1ms, which is much faster than the oscil- 
lation period and slow enough that the BEC remains in 
the lower eigenstate. This process effectively gives the 
BEC a pulsed momentum, named synthetic electric force 
in Ref. [J]. Further, by varying <5 m , one excites the dipole 
oscillation with different amplitudes. The BEC starts to 
oscillate at t\ and is held for a holding time th = £2 — 1\. 
At £2 both the Raman lasers and the trap are quickly 
switched off within 1/is. With the Stern-Gerlach tech- 
nique, a time-of-flight (TOF) image is taken after 24ms 
of free expansion to map out the spin and momentum of 
the BEC. A comparison experiment without SO coupling 
is also carried out by setting a large detuning 8 > 20E I 
instead of turning off the Raman lasers. Accordingly, the 
time sequence for Bbi as is modified: Bbi as is initially set 
for 8 = 8{ at t = and quickly switched to be 8 > 20E r 
at £2 (red line in Fig.QJd)). 

Momentum oscillation. Figure [Ue-f ) displays typical 
momentum oscillations observed in the experiment. The 
frequency along x is significantly changed by the SO cou- 
pling, while those along y and z remain to be the trap 
frequency (independent of oscillation amplitude). Fig- 
ure [2^a) shows oscillation period T x versus Q along path 
Pi (8 — 1E T ). As n increases, the deviation from the 
trap frequency becomes larger. Moreover, it is found 
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that T x is amplitude-dependent, as shown in Fig. [2fb) 
for VL — 4.8EV and 5 = 1E Y . These two features clearly 
demonstrate that the dipole oscillation along the x direc- 
tion is no longer harmonic. 

Similar il-dependence of T x along P^ and P3 are dis- 
played in Fig. [^c) and (d). In the yellow regime of 
Fig-QJc), the single-particle spectrum has two nearly de- 
generate minima separated by a barrier as low as trap- 
ping energy. Macroscopic quantum tunneling between 
two minima takes place in this regime, and the dipole 
oscillation becomes rather complicated and does not fit 
a single-frequency oscillation (shown in the inset of Fig. 
Hd)). 

This anharmonic behavior can be understood from 
the equation-of-motion. Consider Hamiltonian H = 
Zi(Ro,i + (l/2)m^ r 2 () + ^.Uiri - rj), where a 
sums over 2, y and z, U(r) represents two-body inter- 
action. Let X := (l/N)J2i x i be- the center-of-mass 
displacement operator (N be the number of atoms), 

the equation-of-motion for X can be derived as X = 



b 1.0 



1/(N m)Y,i{kx,i+k r a z ,i) and A = -oj 2 x X + iVLk v J2i °~ 
For either no coupling between spins (tt = 0) or 
momentum-independent coupling-e.g. coupled by radio- 
frequency field with k T = 0, the equation-of-motion will 
close at the second order, which yields harmonic oscilla- 
tion and the well-known Kohn theorem [l8|, For both 
nonzero f2 and k r , however, the equation-of-motion can- 
not close at any finite order. This results in anharmonic 
dipole oscillation. 

For a quantitative calculation of the oscillation fre- 
quency, we shall apply a variational wave-function ap- 
proach [26] . We first assume that the condensate stays in 
the lower eigen-branch during the entire oscillation 24 1. 
We further ignore spin-dependent interaction because the 
spin-dependent interaction energy is about 0.46% of the 
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total energy for the F ~ 1 manifold of 87 Rb atoms 
and the aspect ratio of condensate in our experiments is 
far away from "mode resonance" 24j. With these sim- 



plifications, the dipole oscillation is described by [24J 

k x = -u 2 x x , x = d£(k x )/dk x , 



(2) 



As shown in Fig. [5J our experimental data agree well 
with calculations based on Eq. In particular, for 

(il = 4£' r , 5 = 0) where effective-mass approximation 
breaks down because of divergent effective mass, Fig. [2] 
shows that the oscillation frequency remains finite. 

Magnetization Oscillation. The Stern-Gerlach TOF 
images in Fig. (3{a) show that during the dipole oscil- 
lation, the spin population also oscillates. Figure [3Jb) 
and (c) display the oscillation of quasi-momentum k x and 
of polarization M. = (n^ — nj_)/(n^ + n^), respectively; 
it can be seen that their frequency is exactly identical. 
To provide an intuitive picture, we assume that dur- 
ing the entire oscillation, the BEC remains at the lower 
eigenstate branch of Hamiltonian ([1} , and the spin con- 
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FIG. 3: Magnetization oscillation for S — IE T and Q, = 4.8-E r . 

(a) , Spin-resolved TOF images for various holding times th- 

(b) , Quasi-momentum k x versus t^. (c), Polarization M ver- 
sus th- Red circles are directly measured, while black squares 
are deduced from (b) (see text for details). 



figuration adiabatically follows the center-of-mass mo- 
tion. The eigenstate wave function of the lower branch 
"0 = ft(kx)\ t) + fi{k x )\ 4) can be obtained by diagonal- 
izingEq. ©, which yields M(k x ) = \f t (k x )\ 2 '- \ h(k x )\ 2 . 
Hence, the magnetization changes with k x . The function 
M.{t) can be then obtained by using the experimentally 
measured k x (t) (Fig. EJb)). The results are shown as 
black squares in Fig. EJc), and agree very well with the 
directly measured data (red circles). The magnetization 
oscillation reflects the locking between spin and momen- 
tum in Hamiltonian ([1]), and further provides a direct 
justification for the assumption in the variational wave- 
function approach. 

The magnetization oscillation can be understood from 
the absence of Galilean invariance. Consider a BEC that 
moves with velocity v along x, in the co-moving frame the 
single-particle Hamiltonian acquires an additional term 
vk x . In a conventional BEC, this term can be gauged 
away by a gauge transformation ip — >■ e lmvx ip. In our sys- 
tem, however, such a procedure will introduce a velocity- 
dependent Zeeman-energy term — mufc r cr z . Hence, once 
the condensate moves, an oscillation of magnetization Ai 
has to be induced. 

Spin Polarization Susceptibility and Quantum Phase 
Transition. A unique feature of SO coupled condensate 
is that spin polarization susceptibility can be deduced 
from the amplitude ratio between momentum and mag- 
netization oscillations (25[. Here we focus on 5 = 0. For 
n < 4EV, bosons condense in one of the double minima, 
which spontaneously breaks the time-reversal symmetry 
and has non-zero magnetization (o- z ). For D, > 4_E r , 
bosons condense in zero-momentum state with zero mag- 
netization. Thus, there is a phase transition from mag- 
netic condensate to nonmagnetic condensate as f2 varies 
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FIG. 4: Amplitude ratio of spin and momentum oscillation 
Acr/(Ak/kr) versus fl, for magnetic phase (a) and nonmag- 
netic phase (b). The inset is for the spin polarization suscep- 
tibility E-c x deduced from this ratio via Eq. Q . The red 
solid line is from the solution of Eq. ([3]). 



[23l |25| . The spin polarization susceptibility x can be 
expressed as [251 ] 



n 2 /2E, 



2 



for Q < AE r 
for fl > AE t 



(3) 



Equation ([3]) predicts that x diverges at the phase tran- 
sition point D, = AE T . It is further proposed [25| that x 
can be measured from the amplitude ratio between spin 
and momentum oscillations via 



A a = E lX 
A k /k r 1 + E tX 



(4) 



We experimentally measure ratio A a /(A)~/k T ) and com- 
pare the data to the theoretical prediction by Eq. (|3|). 
Figure HJa) and (b) shows the results for magnetic phase 
with ft < 2.5E r and for nonmagnetic phase ft > 4_E r , re- 
spectively. Note that for the reason discussed above, no 
experimental data are available for 2.5E r < < AE r . We 
further deduce susceptibility x y i& Eq. (|4]) from the ex- 
perimental ratio A a /(Ak/k T ), shown in the inset of Fig. 
IH The excellent agreement with the theoretical result 
confirms the unique feature of SO coupled condensate 
that spin susceptibility can be measured via dipole oscil- 
lation. In particular, one finds that as fi — > 4E T + + , x 
does display divergent behavior, giving a strong evidence 
of quantum phase transition. 

In summary, we have experimentally demonstrated 
non-trivial properties of the dipole oscillation for a SO 
coupled BEC. From the experimentally measured dy- 
namics of dipole oscillation, we further display the diver- 
gent behavior of spin polarization susceptibility-a static 
physical quantity-at the quantum phase transition. Be- 
side being a direct experimental observation of unconven- 
tional dynamics in a SO coupled condensate, our quan- 
titative results also provide a bench mark for various re- 
cently developed theoretical approaches. It is expected 
that further study of dynamic behavior would provide a 
powerful tool in probing novel phases of the SO coupled 
BEC, such as stripe superfluid phase 12t 1 1 3 1 - 
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